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Abstract We present an interdisciplinary review of the generalized Cerenkov emis- 
• sion of radiation from uniformly moving sources in the different contexts of classi- 

cal electromagnetism, superfluid hydrodynamics, and classical hydrodynamics. The 
details of each specific physical systems enter our theory via the dispersion law of 
the excitations. A geometrical recipe to obtain the emission patterns in both real and 
wavevector space from the geometrical shape of the dispersion law is discussed and 
applied to a number of cases of current experimental interest. Some consequences 
of these emission processes onto the stability of condensed-matter analogs of grav- 
itational systems are finally illustrated. 



1 Introduction 

The emission of radiation by a uniformly moving source is a widely used paradigm 
in field theories to describe a number of very different effects, from the wake gener- 
ated by a swimming duck on the surface of a quiet lake [1, 2, 3, 4, 5, 6, 7, 8, 9, 10], 
to the Cerenkov emission by a charged particle relativistically moving through a 
dielectric medium [11], to the sound waves emitted by an object traveling across 
a fluid or a superfluid at super-sonic speed [12, 13, 14, 15]. Of course, the radi- 
ated field has a different physical nature in each case, consisting e.g. of gravity or 
capillary waves at the water/air interface, or electromagnetic waves in a dielectric 
medium, or Bogoliubov excitations in the superfluid. In spite of this, the basic qual- 
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itative features of the emission process are very similar in all cases and a unitary 
discussion is possible. 

In the present chapter, we shall present an interdiscipUnary review of this gener- 
alized Cerenkov effect from the various points of view of classical electromagnetism, 
superfluid hydrodynamics, and classical hydrodynamics. The details of each specific 
physical systems enter our theory via the dispersion law £2 (k) of its excitations. In 
particular, the emission patterns in both real and wavevector space can be extracted 
from the geometrical shape of the intersection of the Q (k) dispersion law with the 
Q =k \ hyper-plane that encodes energy-momentum conservation. Once the dis- 
persion law of a generic system is known, our geometrical algorithm provides an 
efficient tool to obtain the most significant qualitative features of wake pattern in a 
straightforward and physically transparent way. 

In the last years, the interest of the scientific community on this classical problem 
of wave theory has been revived by several experiments which have started explor- 
ing the pecuUar features that appear in new configurations made accessible by the 
last technological developments, e.g. the Cerenkov emission of electromagnetic ra- 
diation in resonant media [16, 17] and the Bogoliubov-Cerenkov emission of sound 
waves in bulk dilute superfluids [13, 14, 15]. Another reason for this renewed inter- 
est comes from the condensed-matter models of gravitational systems that are the 
central subject of the present book. In many of such analog models, the presence 
of a horizon may be responsible for the emission of waves from the horizon by 
generalized Cerenkov processes. A full understanding of these classical effects is 
then required if one is to isolate quantum features such as the analogs of Hawking 
radiation, dynamical Casimir emission and anomalous Doppler effect. 

The structure of the chapter is the following. In Sec.2, we shall introduce the 
general field-theoretical formalism to calculate the real and momentum space emis- 
sion patterns and the geometrical construction to obtain qualitative information on 
them. These methods will then be applied in the following sections to a few dif- 
ferent systems of current interest. As a first example, in Sec. 3 we wiU review the 
main features of the Cerenkov emission of electromagnetic waves from relativisti- 
caUy moving charges in a dielectric medium. We shall restrict our attention to the 
simplest case of a non-dispersive dielectric with frequency-independent refractive 
index n, where the phase and group velocities are equal and constant. In this case, 
a Cerenkov emission takes place as soon as the charge speed exceeds the velocity 
of Ught in the medium c = co/n. Modern developments for the case of a strongly 
dispersive media [16, 17, 18], photonic crystals, and left-handed metamaterials [19] 
will be briefly mentioned. In Sec.4 we shall review the emission of sound waves by 
a supersonically moving impurity in the bulk of a dilute superfluid. In addition to the 
Mach cone that appears in the wake behind the object, the presence of single -particle 
excitations in the excitation spectrum is responsible for the appearance of a series of 
curved wavefronts ahead of the impurity. On the other hand, a subsonically moving 
impurity will produce no propagating wave and the perturbation will remain local- 
ized in its vicinity: the resulting frictionless motion is one of the clearest examples 
of the class of phenomena that go under the name of superfluidity [21, 22, 23]. This 
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physics is currently of high experimental relevance, as first real-space images of 
the density perturbation pattern induced by a moving impurity have been recently 
obtained using Bose-Einstein condensates of ultracold atoms [14] and of exciton- 
polaritons in semiconductor microcavities [13]. The case of a paraboUc dispersion 
will be presented in Sec. 5: this specific functional form allows for an elementary 
analytical treatment of the wake pattern in both real and momentum space. On one 
hand, this discussion provides an useful guideline to understand the qualitative shape 
of the wake in superfluids and in surface waves. On the other hand, it is of central 
importance in view of the experimental realization of gravitational analogs based on 
magnon excitations in magnetic soUds [24]. The physics of a material object such 
as a boat, a duck or a fishing line creating surface waves on the air/water interface 
of a lake will be considered in Sec. 6: not only does this example provide the most 
intuitive example of the generalized Cerenkov effect, but is perhaps also the richest 
one in terms of different behaviors that can be observed depending on the system 
parameters, e.g. the velocity of the object with respect to the fluid, the depth of the 
water, the surface tension of the fluid [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]. The concepts that 
have been laid down so far are finally applied in Sec. 7 to analog models of grav- 
ity based on flowing superfluids or surface waves on flowing water. In both these 
cases, classical Cerenkov emission into the so-called zero modes at the horizon may 
disturb detection of the analog Hawking radiation as well as affect the dynamical 
stabihty of the analog black/white hole [25]. Conclusions are finally drawn in Sec. 8. 



2 Generic model 



In this section, we introduce the generic model that will be used to study the dif- 
ferent physical systems in the following sections. The model is based on a Unear 
partial differential equation for a scalar C-number field (p{r,t): in most relevant 
cases, the multi-component physical field (i.e. the vector electromagnetic field or 
the BogoUubov spinor) can in fact be reduced to a single scalar field upon straight- 
forward algebraic manipulations under controlled approximations. We are also as- 
suming that quantum fluctuations of the field <^)(r,f) can be fully neglected. The ge- 
ometry under investigations consists of a spatially homogeneous system interacting 
with a spatially localized moving source describing the moving electric charge, or 
the interaction potential of the moving impurity with the fluid, or the extra pressure 
exerted on the fluid surface by the moving object. In this geometry, the microscopic 
information on the field dynamics is summarized in the dispersion law relating the 
frequency 12 of a plane wave to its wavevector k: different forms of dispersion laws 
corresponding to first- or higher-order partial differential equations are discussed in 
the subsections Sees. 2.1 and 2.3. The geometric construction of the wiike pattern 
starting from the dispersion law £2 (k) is discussed in Sec. 2.2. 
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2.1 The wave equation and the source term 



We start by considering a generic, t/-dimensional classical complex field (p(r,t) (d = 
2 in the figures) that evolves according to the generic linear, first-order in time, 
partial differential equation: 

!a,0(r,f) =i2(-iVr)0(r,O+5(r,O (1) 

with a source term S{r,t). 

The function Q (k) defines the so-called dispersion law for free field propagation, 
that is the frequency of the plane wave solutions 

<j){r,t) = <l)oe''^e-'^^^^'. (2) 

as a function of the wavevector k in the absence of sources, S{r,t) = 0. 

Throughout this chapter we shall consider a spatially localized and uniformly 
moving source term of the form 

S{r,t) = So{r-yt), (3) 

with a spatial profile 5'o(r) concentrated in the vicinity of r = and moving at a 
speed v. 

Thanks to the translational invariance of the free field problem in both space and 
time, solution of the full wave equation (1) in the presence of the source term is 
easily obtained in Fourier space with respect to both space and time. Defining the 
Fourier transform in the usual way 

^ (k, £») = jdtjd'^r(j>{r,t) g-*-" e''" , (4) 

the source term in Fourier space has the simple form 

5(k,ffl) =2;r5o(k)5(»-k-v) (5) 

in terms of the structure factor >So(k) defined as the Fourier transform of the source 
shape 5o(r). 

In Fourier space, the solution of (1) is then 

, 27tSo(k)S{co-k-y) 
'^('^'^^= co-X2(k) + ,-0+ ' 

where an infinitesimally small imaginary part is introduced in the denominator of 
(6) to specify the integration contour to be followed around the poles and, in this 
way, ensure causality of the solution. This trick dates back to Rayleigh [4] and is 
equivalent to a infinitesimal shift of the dispersion law into the lower half-space, 
X2 (k) — ^ X2 (k) — iO^. Physically, it corresponds to introducing a very weak damping 
of the plane wave solutions in time. 
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Hr,t)=(^oe'^e-'''^''^'e-'^*'. (7) 

or to assume that the source term is slowly switched on in time [9]. 

The real-space pattern is obtained by an inverse Fourier transform of (6), 

Thanks to the 5(0) — k • v) factor in (6), this expression only depends on the com- 
bination r' = r — V/: as expected, the wake pattern is rigidly moving at the speed of 
the source. Within GaUlean invariance, the r' coordinate corresponds to the spatial 
coordinate in the reference frame of the source in motion at velocity v. 

Evaluation of (8) can be performed with standard numerical tools. The result 
for some most interesting cases will be presented in the next sections. Now, we 
shall rather proceed with some analytical manipulations of (8) that allow to extract 
qualitative information on the emitted field pattern from the dispersion law 12 (k). 
The first step in this direction is to note that the integral in (8) is dominated by those 
k values for which the resonant denominator vanishes, that is 

i2(k)=k-v. (9) 

This equation recovers the standard Cerenkov condition for emission of radia- 
tion [11] and can be geometrically interpreted as the intersection of the dispersion 
surface Q (k) with the i2 = k • v plane. In a quantum description of the Cerenkov 
emission by a massive charged particle, the condition (9) naturally appears when 
energy-momentum conservation is imposed to the photon emission process [11]. 
When reformulated in the reference frame of the moving source, the condition (9) 
reduces to Q' — y{Q — k • v) =0, meaning that the perturbation pattern around the 
source at rest is stationary in time in the moving reference frame. 

The locus of k ^ modes that satisfy (9) is a central object in all the following 
discussion as it defines the modes in k space into which the Cerenkov emission 
will be peaked ^ . In particular, no emission of propagating waves takes place if the 
locus Z is empty; the non-resonant contributions to (8) only provide a non-radiative 
perturbation that remains spatially localized in the close vicinity of the source and 
is not able to transport energy away. In spite of this, the momentum and energy that 
are stored in the localized moving pattern of the field (j) are responsible for a sizable 
renormahzation of the particle mass [9, 27, 28]. 



' The k = mode corresponds to a spatially constant modulation that does not transport energy 
nor momentum. As discussed in [26], many other interesting features of wave propagation can be 
graphically studied starting from isofrequency surfaces analogous to the locus E. 
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2.2 Qualitative geometrical study of the wake pattern 

Let us consider a generic point ko E E. Within a neighborhood of ko, we introduce 
a new set of k-space coordinates defined as follows: for each point q, q„ is the 
distance of k from the L surface and the position of the closest point on the surface 
Z is parametrized by the {d — 1) -dimensional q|| curvilinear coordinate system. A 
sketch of this coordinate system is indicated as a grid in Fig. 1(a). 




Fig. 1 Left (a) panel: k-space sketch of a patch of the locus L around the wavevector ko. The 
grid shows the (qy , q„) coordinate system used in the geometrical construction of the wake pattern. 
Right (b) panel: sketch of the region of the wake pattern generated by the emission in the neigh- 
borhood of ko: the blue fringes have wavevector ko, the direction of propagation is determined 
by the normal to the locus E at the point ko. 



In this new coordinate system, the Fourier integral giving the emitted field pattern 
in real space can be approximately rewritten as: 



where k(q|| ) is the position of the point on the surface E corresponding to the value 
q|| of the {d ~ 1) -dimensional coordinate and n is the unit vector normal to the 
surface at kg in the direction of growing 12 (k) — k • v. As the surface E is defined by 
the zeros of i2 (k) — k • v, the velocity 

v; = v/n = Vk[f2(k)-k-v] =vg-v (11) 

is directed along the normal n and corresponds to the group velocity of the wave, 
as measured relative to the moving source at v. For a non-relativistic source speed 
V <C Co, it can be interpreted as the group velocity observed from the source reference 
frame. 

The integral over qn can be performed by closing the contour on the complex 
plane. The only pole is located slightly above the real axis. Depending on the sign 
of n • r', the contour has to be closed in the upper or lower half plane, which gives 
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^{r') = -iS,{k,) j^-^,^^@W,-A- (12) 

The expression (12) can be further simplified by performing the so-called sta- 
tionary phase approximation, as first proposed by Thomson [4]. For each value r' of 
the relative coordinate, the integral over q|| is dominated by those points for which 
the phase is stationary, i.e. the variation of k(q|| ) on q|| is orthogonal to r'. In com- 
bination with the Heaviside-0 function in (12), this is equivalent to requiring that 
the vector r' is parallel to the normal n to the surface X, at point ko in the direction 
of growing i2 (k) — k • v , i.e. parallel to the relative group velocity v^. 

For a generic relative position f ', there are only a few discrete points k^ on L such 
that this condition is met. As a consequence, for generic values of r' in a neighbor- 
hood of r', one can approximately write 

where the sum is over the allowed k, vectors: the numerical coefficient Ak^j is in- 
versely proportional to the curvature of I, at k^ and v|, j is the group velocity of the 
kj mode. 

A physical understanding of this result can be easily obtained by looking at 
the diagram of Fig. 1(b). Within Galilean invariance, sitting in the moving refer- 
ence frame of the source may facihtate building an intuitive picture of the emission 
process: every point on the surface Z corresponds to a continuous plane wave of 
wavevector ko that is emitted from the source and propagates away from it at a 
group velocity v^, (indicated by the green arrow in the figure). As a result, it is able 
to reach all points r' that lie in the vicinity of the straight line of direction v^,. While 
the group velocity is always along the radial direction, the wavevector ko (blue ar- 
row in the figure) can have arbitrary direction: as a result, the wavefronts (indicated 
by the blue fringes) are generally tilted and the emission pattern does not necessarily 
resemble to a spherical wave. Of course, all this reasoning can be performed equally 
well in the laboratory frame if is interpreted as the relative group velocity of the 
wave with respect to the moving source. 



2.3 Generalization to higher-order wave equations 

The geometrical framework introduced in the previous subsections is not restricted 
to partial differential equations that are of first-order in time, but can be extended to 
more general wave equations of the form 



P[idt-Nr]^{Y,t)=S{Y,t), 



(14) 
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where P is an arbitrary polynomial in two variables, a scalar variable and a d- 
component vectorial variable. The degree of the polynomial P corresponds to the 
order of the partial differential equation for <^>(r,/): in the case of electromagnetic 
waves in a non-dispersive medium, it is of second order in both variables; in the case 
of Bogoliubov excitations in a superfluid, it is of second order in time and of fourth 
order in space; in the case of surface waves, it is of second order in time, but it in- 
volves arbitrarily high derivatives in the spatial coordinates. The different branches 
Q, (k) of the dispersion law are then defined by the roots of P via the equation 

P[X2(k),k]=0. (15) 

In the presence of a source term of the form (3), the solution of (14) has the form 

where the infinitesimally small imaginary part has been again added in order to 
enforce causaUty by shifting the real roots Q, of the dispersion law (15) into the 
lower half of the complex -plane. 

The reasoning to extract from (16) the qualitative shape of the real-space pattern 
is then the same as before, the locus X in k-space being now defined by the zeros of 
the polynomial equation 

P(k-v,k)=0 (17) 

with k 7^ 0. In the next sections, we shall discuss in full detail a few physical ex- 
amples illustrating how the geometrical structure of T, determines the shape of the 
emission pattern in both real and momentum spaces. 



3 Cerenkov emission by uniformly moving charges 

As a first application of the theory, in this section we shall review the basic features 
of the emission of electromagnetic radiation by a charged particle relativistically 
moving through a dielectric medium at speed higher than the phase velocity of light 
in the medium. This is the well-known Cerenkov effect (or, more precisely, Vavilov- 
Cerenkov effect) first observed by Marie Curie, then experimentally characterized 
by Vavilov and Cerenkov [29] and finally theoretically understood by Frank and 
Tamm [30]. 



3.1 Non-dispersive dielectric 



In the simplest case of a non-dispersive dielectric with a frequency-independent 
refractive index n, the dispersion law satisfies the second-order equation 
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(18) 



in the (^2 , k) space, the dispersion Q (k) corresponds to a conical surface with vertex 
in 12 = A: = 0. A cut of this cone along the ky = Q line is shown in Fig. 2(a): the 
thick lines indicate the positive frequency part of the conical surface, the thin line 
indicate the negative frequency part. In the absence of dispersion, the phase and 
group velocities coincide and are equal to c = co/n. 

The shape of the locus Z of k ^ points satisfying i2 (k) — k • v = crucially 
depends on whether the charge is moving at a sub-luminal v < c or super-luminal 
V > c speed. In the former case, the locus L is empty and no radiation is emitted. 
The localized, non-radiative perturbation that is present around the charge due to the 
non-resonant excitation of the field modes contributes to the (velocity-dependent) 
Coulomb field around the charge. 

The locus i; in the case of a super-luminally moving charge in the positive x 
direction is illustrated in Fig. 2(b): it has the analytic form 



and consists of a pair of half straight lines originating from k = and symmetrically 
located with respect to the kx axis at an angle Q such that cos Q = c/v. The higher 
the particle speed v/c, the wider the angle d made by the direction of the Cerenkov 
emission with the direction of the charge motion. 

The most pecuUar feature of the locus I. is that the normal vector to Z indicating 
the direction of the relative group velocity v^ = [Q, (k) — k • v] = v^ — v [indicated 
by the green arrows on Fig.2(b)] is constant for all points k lying on each of the 
two straight lines forming Z and points in the backward direction. This last feature 
is a straightforward consequence of the fact that the charge velocity is larger than 
the speed of light c in the medium. As a result, all modes on Z propagate (as seen 
from the charge reference frame) into the same direction and the electromagnetic 
field radiated by the charge is spatially concentrated around the direction of v^. This 
defines a single-sheet conical surface in real space, i.e. a pair of half straight lines in 
the two dimensional geometry considered here. 



Its aperture around the negative x axis ^ is determined by the condition sin = c/v: 
the faster the charge speed, the narrower the cone behind the charge. In the analogy 
with the conical sonic wake generated by a super- sonically moving bullet in a bulk 
fluid, we will refer to this real space cone as the Mach cone. The very thin shape 
of the Mach cone results from the interference of the continuum of points on the 
E locus. For each ko € E, the fringes are orthogonal to the Mach cone and have 




(19) 



2 C X 

y = 2 with x<0 



(20) 



^ The coefficients of the analytical form (20) can be understood from the Fourier ttansform of a 
delta function peaked on the conically-shaped locus E of equation (19). 
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Fig. 2 Panel (a): cut along ky = of the photon dispersion in a non-dispersive medium of 
frequency-independent refractive index n. The dashed line indicates the Q. = k-v plane for a super- 
luminal charge speed v > c. Panel (b): k-space locus L of resonant modes into which the Cerenkov 
emission occurs, the so-called Cerenkov cone; the green arrows indicate the nornial to the Z locus, 
that is the direction of the relative group velocity vj, = \fi (k) — k • v] . Panel (c): real-space pat- 
tern of the electric field amplitude in the wake of the charge; the pattern is numerically obtained as 
the fast Fourier transform of the k-space perturbation (16). The Mach cone around the negative x 
axis is apparent with aperture . 
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different spacing: the interference is everywhere destructive but for the thin surface 
of the Mach cone. If the correct form of the structure factor >S'o(k) is included, the 
usual 5-shape for the Mach cone is recovered [11]. 

In view of the discussion of the next sections, it is crucial to clearly keep in mind 
the conceptual distinction between the Mach cone in real space on which the elec- 
tric field intensity is spatially concentrated and the k-space Cerenkov cone defining 
the directions into which the radiation does occur. The former was experimentally 
detected and characterized in [17, 31] by looking at the spatial profile of the electric 
field wake behind the charge ^. The latter is observed in any standard Cerenkov ra- 
diation experiment measuring the far-field angular distribution of the radiation, that 
turns out to be concentrated in the forward direction on a conical surface making a 
Cerenkov radiation angle 9 with the charge velocity. 

The conceptual distinction between the Cerenkov and the Mach cones is related 
to the distinction between the so-called phase and group cones, first pointed out in 
the context of the Cerenkov emission in dispersive media in [32]. Restricting for a 
moment our attention to a given emission frequency, the wave cone is defined as the 
real space conical wavefront passing through the source and orthogonal to the direc- 
tion of the far-field emission in k space: its aperture <^p/j around the negative x axis is 
determined by the phase velocity as sin t^p^ = Vph/v and is related to the aperture of 
the Cerenkov cone by ^p/, = 71/2—6. With some caveats, it can be interpreted as the 
wavefront on which the Cerenkov emission has a constant phase. On the other hand, 
the group cone is defined as the Mach cone for the given frequency and describes 
the spatial points on which the (spectrally filtered) electric field intensity is peaked. 
Its aperture ^ depends on the group velocity of light Vgr as sin^ = Vgr/v. The dis- 
tinction between the phase and group cones has been anticipated in [18] to be most 
striking in the case of ultra-slow light media where Vgr is reduced to the m/s range 
while Vph remains of the order of the speed of light in vacuo cq — 3 • 10^ m/s [33]. 

The study of the Cerenkov effect in strongly dispersive media where the refrac- 
tive index n{03) has a strong dependence on the frequency and/or the medium ex- 
hibits a non-trivial spatial patterning is still a very active domain of research from 
both the theoretical and the experimental points of view. For instance, the conse- 
quences of a strong resonance in n{(o) were theoretically investigated in [16]: the 
sub-linear dispersion of the photon in a resonant medium is responsible for the dis- 
appearance of the threshold velocity for the Cerenkov emission and for a non-trivial 
spatial patterning of the electric field wake behind the charge. These striking results 
were experimentally confirmed in [17] and bear a close resemblance to the surface 
wave pattern in the wake of a duck swimming on shallow water that will be dis- 
cussed in Sec. 6. 3. Another active and promising research line is addressing those 
new features of Cerenkov radiation that follow from the peculiar band dispersion of 
photons in spatially periodic media [19] and in negative refractive index metamate- 
rials, the so-called left-handed media [20]. 

' It is interesting to note that in both these experiments the moving charge responsible for the 
Cerenkov emission did not consist of a charged physical particle travelling through the medium, 
but rather consisted of a moving bullet of nonlinear optical polarization generated by a femtosecond 
optical pulse via the so-called inverse electto-optic effect. 
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4 Moving impurities in a superfluid 



A central concept in the theory of superfluids [21, 22, 23] is the so-called Landau 
criterion for superfluidity, that determines the maximum speed at which a weak 
impurity can freely travel across a superfluid without experiencing any friction force 
and without generating any propagating perturbation in the fluid. In terms of the 
dispersion Q (k) of the excitations in the superfluid, the Landau critical velocity has 
the form 



mm 
k 



(21) 



This cornerstone of our theoretical understanding of quantum liquids has a simple 
interpretation in terms of the theory of the generalized Cerenkov effect reviewed in 
Sec. 2: the friction force experienced by the moving impurity is due to the emission 
of elementary excitations in the fluid by a mechanism that is a quantum fluid ana- 
log of Cerenkov emission. The v < Va- condition for superfluidity corresponds to 
imposing that the locus E of excited modes is empty, while for faster impurities a 
characteristic wake pattern is generated around the impurity. 




Fig. 3 Left panel: experimental image of the real-space wake pattern that appears in a Bose- 
Einstein condensate of ultracold atoms hitting the repulsive potential of a blue-detuned laser beam. 
The condensate motion is from right to left. Figure taken from [14], as published in [15]. Middle 
and right panels: experimental images of the real-space wake pattern (middle) and the momentum 
distribution (right) for a Bose-Einstein condensate of exciton-polaritons hitting a fabrication defect 
in the planar microcavity. The polariton flow is from top to bottom. The value of the density in the 
right panel is very small and interactions negligible. Figures taken from [13]. 



An experimental image of this wake using a dilute Bose-Einstein condensate of 
ultracold atoms hitting^ the repulsive potential of a blue-detuned laser is reproduced 
in the left panel of Fig. 3; an analogous image for a condensate of exciton-polaritons 
in a semiconductor microcavity is reproduced in the middle panel. In both cases, the 
density wake extends both behind and ahead of the impurity. The geometrical shape 

* Needless to say that the configuration of a moving superfluid hitting an impurity at rest is fully 
equivalent modulo a Galilean transfomiation to the case of a moving impurity crossing a superfluid 
at rest. 
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of the Z locus is instead clearly visible in the momentum distribution pattern shown 
in the right panel. 

The situation is of course more complex when stronger impurities are considered, 
e.g. a finite-sized impenetrable object: in this case, the critical speed for frictionless 
flow was predicted in [34] to be limited by the nucleation of pairs of quantized 
vortices at the surface of the object, and therefore to be significantly lower than 
the speed of sound. This mechanism was recently confirmed in experiments for 
with atomic [35] andpolariton [36] condensates. Furthermore, it is worth reminding 
that all our reasonings are based on a mean-field description of the condensate that 
neglects quantum fluctuations: more sophisticated Bethe ansatz calculations for a 
strongly interacting one-dimensional Bose gas [28] have anticipated the appearance 
of a finite drag force also at sub-sonic speed. Including higher order terms of the 
Bogoliubov theory led the authors of [37] to a similar claim for a three-dimensional 
condensate. 



4.1 The Bogoliubov dispersion of excitations 

The theoretical description of superfluids is simplest in the case of a dilute Bose gas 
below the transition temperature Tbec for Bose-Einstein condensation [21]. For T <C 
Tbbc and sufficiently weak interactions, most of the atoms are accumulated in the 
same one-particle orbital, the so-called Bose-Einstein condensate. The elementary 
excitations in a dilute Bose-Einstein condensate are characterized by the Bogoliubov 
dispersion [21] 

h^Q^ = — [—+2^], (22) 

where m is the mass of the constituent particles and the chemical potential is given 
(at zero temperature) by 

M = (23) 

m 

in terms of the particle-particle low-energy collisional scattering length ao and the 
particle density n. In the standard three-dimensional case, the weak interaction (or 
diluteness) condition requires that hoq ^ 1 

The characteristic shape of the Bogoliubov dispersion (22) is illustrated in the 
(a,d) panels of Fig.4. For small momenta k^^ (the so-called healing length ^ 
being defined as /i'^/wi^^ = ;U), the dispersion has a sonic behavior 

X22~c?A^ (24) 

with a soimd speed Cs = s/jijm, while at high wavevectors ^ I it grows at a 
super-sonic rate and eventually recovers the paraboUc behavior of single particles, 

■ (25) 

2m 
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Fig. 4 Top row: Bogoliubov dispersion of excitations in a dilute Bose-Einstein condensate. The 
dashed line indicates the X2 = k • v plane for two different impurity speeds v/cj = 1.2 [panel (a)] 
and v/Cs = 2.5 [panel (d)]. Middle row: shape of the corresponding k-space locus Z of resonantly 
excited modes. The dashed lines indicate the Cerenkov cone in the low wavevector region ki, <Si\; 
the green an'ows indicate the nornial to the L locus, that is the direction of the relative group 
velocity v^. Bottom row: real space pattern of the density modulation. All patterns are numerically 
obtained performing the integral via a fast Fourier transform of the k-space perturbation (16). The 
black dashed lines indicate the Mach cone. 
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An explicit calculation from (22) shows that the Landau critical velocity (21) in the 
dilute Bose gas is determined by the speed of sound v^r = c^. It is worth remind- 
ing that this is no longer true in more complex superfluids with strong interparticle 
interactions as liquid He-II, where Vcr is determined by the roton branch of the el- 
ementary excitations [21, 22, 23, 38]. Remarkably, super-linear dispersions in the 
form (22) also appear in the theory of surface waves on shallow fluids when the 
fluid depth is lower than the capillary length, see Eqs.42 in Sec.6. 

The effect of the moving impurity onto the superfluid can be described by a time- 
dependent external potential of the form V{r,t) =Vo{r — \t) coupled to the particles 
forming the superfluid. Inclusion of this external potential in the Bogoliubov theory 
requires including a classical source term in the Bogoliubov equations of motion 
for the two-component spinor describing the quantum field of the non-condensed 
particles: a complete theoretical discussion along these Unes can be found in the 
recent works [12, 15]. Here we shall use an approximate, yet qualitatively accurate 
model based on the simplified scalar theory of Sec. 2: the real and imaginary parts of 
the field ^{r,t) correspond to the density and phase modulation of the condensate. 



4.2 Superfluidity vs. Bogoliubov-Cerenkov wake 

As we have already mentioned, the locus Z is empty for sub-sonic impurity speeds 
V < Cs'. the impurity is able to cross the superfluid without resonantly exciting any 
propagating mode of the fluid. As a result, within mean-field theory it is not expected 
to experience any friction force. Modulo a Galilean transformation, this effect is 
equivalent to a frictionless flow along a containing pipe in spite of the roughness of 
the walls, which is one of the clearest signatures of superfluid behavior [21, 22, 23]. 

Still, the non-resonant excitation of the BogoUubov modes by the moving impu- 
rity is responsible for a sizable density modulation in the vicinity of the impurity, 
that quickly decays to zero in space with an exponential law. An important conse- 
quence of this localized density perturbation is a sizable renormalization of the mass 
of the object [28]: the linear momentum that is associated to the moving impurity 
gets in fact a contribution from the portion of fluid that is displaced by it. 

For super-sonic motion, the locus Z consists of a conical region at small 
analogous to what was found in Sec. 3.1 for a purely linear dispersion: as in that 
case, the aperture angle 9 of the k-space Cerenkov cone [dashed lines in Fig.4(b,e)] 
defining the far-field angle at which phonons are emitted by the impurity is defined 
by the condition cos = c^j v. 

Correspondingly, the aperture ^ of the Mach cone that is visible in the real-space 
density modulation pattern behind the impurity is defined by sin0 = Cj/v (dashed 
black lines in [Fig.4(c,f)]. This cone is the superfluid analog of the Mach cone that 
is created in a generic fluid by a super- sonically moving object, e.g. an aircraft or a 
bullet. An experimental image of a Mach cone in a superfluid of exciton-polaritons 
is shown in the central panel of Fig.3. As usual, the faster the impurity, the narrower 
the Mach cone. 
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Differently from sound waves in an ordinary fluid, the Bogoliubov dispersion of 
the excitations in a superfluid is characterized by a paraboHc shape at large wavevec- 
tors ^ I according to (25). This region of the BogoUubov spectrum is respon- 
sible for the smooth arc in the high wavevector region of Z that connects the two 
straight lines emerging from the origin k = 0. In experiments, the shape of E can 
be inferred following the peak of the momentum distribution of the particles in 
the superfluid: an example of experimental image using exciton-polaritons in the 
low-density regime is reproduced in the right panel of Fig.3. Analogous images for 
atomic gases can be found e.g. in [39]. 

In the low wavevector region, the relative group velocity is oriented along the 
edges of the Mach cone. Along the high wavevector part of E, the relative group 
velocity v^ rotates in a continuous and monotonous way spanning all intermediate 
directions external to the Mach cone. As no point on E corresponds to a relative 
group velocity oriented in the backward direction inside the Mach cone, the density 
profile remains unperturbed in this region. On the other hand, the density pertur- 
bation shows pecuUar features in front of the Mach cone, with a series of curved 
wavefronts extending all the way ahead of the impurity. These wavefronts are clearly 
visible in the experimental images that are shown in Fig.3 for atomic (left panel) and 
polaritonic (middle panel) superfluids. In the k-space diagrams of Fig.4(b,e), these 
waves correspond to the regions in the vicinity of the extreme points of E where 
is directed in the direction of the impurity speed along the positive x direction. 

Physically, these curved wavefronts in the density modulation pattern can be un- 
derstood as originating from the interference of the macroscopic coherent wave as- 
sociated to the Bose-Einstein condensate with the atoms that are coherently scat- 
tered by the moving impurity. An analytic discussion of their shape is discussed in 
detail in [41]; their one-dimensional restriction was first mentioned in [40]. In the 
next section we shall present analytical formulas for an approximated theory where 
the single particle region of the Bogoliubov dispersion is modeled with the parabolic 
dispersion of single-particle excitations. 



5 Parabolic dispersion: conies in the wake 

Another example of dispersion that is fuUy amenable to analytic treatment is the 
parabohc one, 

n{k) = —+n. (26) 
2m 

In spite of its simplicity, this form of dispersion can be used to model a number of 
different physical configurations, from the large wavevector ^ I region of the 
Bogoliubov dispersion (22) of superfluids, to the resonant Rayleigh scattering in 
planar microcavities [42], to magnons in solid-state materials [24, 43]. In particular, 
the results of this section wiU shine further light on the curved wavefronts observed 
in Fig.4(c,f) ahead of the impurity. 




Fig. 5 Top row: parabolic dispersion of excitations (O = /2 + fl in the n = (a), = — 1 < (b), 
fi = 0.4 > (c) cases (for notational simplicity, we have set in = ti= 1). The dashed line indicates 
the i2 = k - V plane for a generic particle speed v = 1 along the positive x direction. Middle row: 
circular shape of the k-space locus X of resonantly excited modes. Bottom row: real space patterns 
of the density modulation. All patterns are numerically obtained performing the integral in (8) via 
a fast Fourier transform algorithm. 



For a generic dispersion of the parabolic form (26), simple analytical manipula- 
tions show that the locus Z has a circular shape as shown in Fig.5(b,e,h). Assuming 
again that the particle speed v is directed along the positive x axis, the center of the 
circle is located at kx = ko = mv/fi, ky — and has a radius k such that 

hi? mv^ 
2m 2h ^ 

Depending on the relative value of the velocity v and of the ji parameter, different 
regimes can be identified. 

For positive (but such that the RHS of (27) is still positive), the radius k is 
smaller than ko and the origin k = lies outside the circle. This is the typical case 
of large k excitations in superfluids, whose dispersion is approximated by Eq.25. 
The usual resonant Rayleigh scattering ring [42] passing through the origin k = 
is recovered in the jj. = case describing the case of an ideal gas of non-interacting 
particles: an experimental example of such a ring is visible in the momentum distri- 
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bution shown in Fig. 3(c) for a low-density gas of (almost) non-interacting polaritons 
flowing against a localized impurity potential. For negative ^, the radius is instead 
larger k> ko and the origin k = falls inside the circle. 

The relative group velocity is directed in the outward radial direction. As a 
consequence of the smooth shape of E, v'^ spans all possible directions and the real- 
space perturbation shown in Fig.5(c,f,i) extends to the whole plane. However, the 
wavefronts can have very different shapes depending on the relative value of v and 

A closed form for the real-space wake pattern is straightforwardly obtained by 
noting that the integral in the right-hand side of (8) is in this case closely related to 
the retarded Green's function for a free non-relativistic particle [44], 



'^'■'''^ = 7 (2^ftkV2m-a)-iO+- ^^^^ 



In a generic dimension d, the asymptotic form of Gret at large r has the outgoing 
spherical wave form 

Gr.M = j^§2^''' (29) 
with a wavevector k such that 

^ = (0. (30) 
2m 

Of course, for co > (or © < 0), the solution such that ^ > (or Im[^] > 0) must be 
considered. Q is a dimension- and energy-dependent normaUzation constant. 

Using this result, the expression (8) for the wake generated by a point-like source 
term can be simpUfied into 

d^k e'^"^ 



A( l\ f ^ 



(27^)2 M + ^-k-v-!0+ 

- y(2.)^^_(«._,)_,o. V? ■ ' ' 

The real part of this wave provides the wake pattern plotted in Fig.5(c,f,i): The 
different panels correspond to the /i = (c), jU < (f) and ;U > (i) cases, which 
correspond tok = k„ (c), k > k„ (f), and k < k„ (i), respectively. 

The shape of the wavefronts is obtained as the constant phase loci of (31). For 
instance, the loci of points for which the phase of the field ^ equals an integer 
multiple of 2n are described by 

k\/x^+y'^ +koX = InM (32) 

with M a generic integer. After moving the kox term to the LHS and then tak- 
ing the square of both members, this equation is straightforwardly rewritten as a 
quadratic equation in the spatial coordinates. The shape of the wavefronts in the 
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two-dimensional plane is therefore described by conic curves: the specific nature of 
the conic in the different cases depends on the ratio k/k,,. 

For kg = k, the wavefronts have a parabolic shape described by the equation 

An^M^ - AnMkoX = koj^ . (33) 

As the square root has by definition a non-negative value, the further condition 
2kM — kgX > has to be imposed to ensure that the RHS of (32) is non-negative. 
Combined with (33), this condition is equivalent to imposing that the integer M > 0. 
An example of these parabolic wavefronts is shown in Fig.5(c). 

For k> ko, the wavefronts have an elliptic shape described by the equation 



-ky+{k'-kt) 



InMko 

X ' 



4K^M^k'^ 

'''' 



An example of these elliptic wavefronts is shown in Fig. 5(f). The condition on the 
non-negativity of the RHS of (32) imposes that M > 0. The ellipticity of £ is a 
function of the ratio k/ko'. the closer this ratio is to 1, the more elongated the ellipse 
is. In the limit k/k,, 1, the ellipse tends to a parabola, recovering the case k = k,, 
discussed above. The larger the ratio k/kg, the closer the wavefront shape to a series 
of concentric circles. 

Finally, for k<ko, the wavefronts have hyperbolic shapes described by the equa- 
tion 

27tMko 



{ki-k') 



-2 2 4;r2M2F 



In this case, the condition on the RHS of (32) does not impose any condition on M 
that can have arbitrary positive or negative integer values. However, combining this 
condition with the equation defining the hyperbola, one finds that for each M only 
the left branch of the hyperbola at lower x has to be retained. Positive vs. negative 
values of M are responsible for the different periodicities that are visible in Fig.5(i) 
in the x > and x<0 regions, respectively. 

From the k space diagrams in Fig.5(b,e,h), it is immediate to see that the points 
on E corresponding to the backward and forward propagating waves are the two 
intersections of the circle with the x axis: the wider spacing of the wavefronts in 
the backward direction is due to the smaller magnitude of the wavevector at the 
intersection point at lower x. In the parabolic case, this point coincides with the 
origin, which explains the absence of density oscillations on the negative x axis. 

The characteristic curvature of the forward propagating wavefronts provides a 
qualitative explanation for the shape of the density modulation experimentally ob- 
served ahead of the impurity and illustrated in the left and central panel of Fig. 3. 
Of course, the absence of backward propagating waves in the superfluid behind the 
impurity is due to the k = singularity of the E locus for the case of the BogoU- 
ubov dispersion. It is worth reminding that, in contrast to previous works, the shape 
of the forward propagating wavefronts is exactly parabolic only in the = limit 
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of non-interacting particles. For the generic n >0 case of Bogoliubov theory, the 
Hartree potential in (25) makes their shape to be closer to (part of) an hyperbola. 



6 Surface waves on a liquid 

The discussion of the previous sections on the Cerenkov effect in classical electro- 
magnetism and on the response of superfluids to moving impurities puts us in the po- 
sition of getting an easy quaUtative understanding of the surface waves that are gen- 
erated by a duck steadily swimming on the surface of a quiet pond or, equivalently, a 
fishing fine in a uniformly flowing river ^. This system is by far the most accessible 
from the experimental point of view, but perhaps also the richest one for the variety 
of different behaviors that can be observed depending on the system parameters. A 
few examples of experimental pictures are shown in Fig. 6. For the sake of concise- 
ness, we shall restrict ourselves to the case of the water-air interface and restrict to 
the linear regime of wave propagation described by the model equation (1). More 
complete treatments based on the full hydrodynamic equations including nonlinear 
effects can be found in the dedicated Uterature, see e.g. [1, 2, 3, 4, 6, 7, 8, 5, 9, 10]. 



6.1 Dispersion of surface waves 



The dispersion of surface waves on top of a fluid layer of height h and at rest has the 
form 

Q {kf = ^gk + ^k^^ tanhkh, (36) 

where p is the mass density of the fluid, g is the gravitational acceleration, and 7 is 
the surface tension of the fluid-air interface. 

In the simplest case of a deep fluid, the tanh A:/? factor can be approximated with 
1 and the dispersion is characterized by two regions. For low wavevectors k <C ky, 
the dispersion follows the sub-hnear square-root behavior 

i2(A:)~±V^, (37) 

of gravity waves, while for large k^ ky it is dominated by capiUarity effects and 
has a super-linear growth as 

n{k)c:^±^k^/^. (38) 



^ It is interesting to note tliat, as in tlie case of electromagnetic waves, accelerated objects emit 
surface waves independently from their speed [48]. 
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Fig. 6 Upper panels: picture of the Kelvin's ship-wave pattern behind a duck swimming at uni- 
form speed on a quiet lake (left). Photograph courtesy of Fabrice Neyret, ARTIS-CNRS, France. 
Experimental picture of the Mach cone downstream of a wire immersed in radially flowing silicon 
oil (upper right). Capillary waves are not visible as they are quickly damped by the larger viscos- 
ity of silicon oil. Picture from [45]. Experimental picture of the Mach cone downstream of a pin 
immersed in very shallow flowing water: the surface wave dispersion is supersonic and the height 
modulation stays outside the Mach cone. Picture courtesy of Silke Weinfurtner (lower right). Mid- 
dle panel: Original hand drawing by John Scott Russell [46] of the waves generated by a vertical 
rod (diameter =1/16 inch) moving along the water surface with a uniform velocity. The rod moves 
in the leftward direction: the capillary waves are visible in front of the rod and the gravity waves 
behind it. A cut of the surface height modulation is shown right above the main drawing. Lower 
panel: Original hand drawing by Lord Kelvin of Kelvin's ship-wave pattern [47]. The BCD wave- 
front belongs to the so-called transverse wave pattern. The so-called diverging waves connect the 
object at A to the dashed lines indicating the edges of the pattern. 



22 



lacopo Carusotto and Gemain Rousseaux 



The characteristic wavevector scale separating the two regions is fixed by the capil- 
lary wavevector 

ky = , 1^. (39) 

V y 

For the specific case of water/air interface, ky ~ 370 m \ which corresponds to the 
value 

iy = I /ky = 2.7 -IQ-^m. (40) 

for the capillary length. 

In fluids of finite depth, one can no longer approximate the tanh in (36) with 1. 
As a result, the dispersion in the low-wavevector region recovers a sonic behavior 
at low k's 

nik)c^±Csk (41) 

with a speed of sound Cg = \/gh proportional to the square root of the fluid depth. 
The sign of the first correction to the sonic behavior (41), 



Q.{kf:^ghk^ 



3 



elk" (42) 



critically depends on the depth of the fluid as compared to the capillary length (40). 
For relatively deep fluids such that h > V^ij, the dispersion has a sub-hnear behav- 
ior, while it recovers a super-linear behavior analogous to the Bogohubov dispersion 
(22) for very shallow fluids such that h < \/lity. 

Independently of the fluid depth h, the super-hnear behavior of the dispersion 
Q. (k) at large k makes the locus Z to be either empty or to consist of a closed 
curve. For very slow sub-sonic motions v < Vmin (with Vmin to be defined in the next 
subsection), the locus Z is empty. For intermediate Cj > v > (or infinitely deep 
fluids, Cj — °°), the locus Z shown in Fig.7(e) consists of a smooth closed curve that 
does not encircle the origin point k = 0. For super-sonic motions v > Cj, the locus 
Z shown in Fig.8(b,e) develops a conical Cerenkov singularity at k = 0. 



6.2 Deep fluid 

Let us start by investigating the deep fluid regime /i — >■ oo for which the sonic speed 
Cs °°. The structure of the locus Z can be understood by looking at Fig.7(d); for 
low speeds 

V<V„,in= (^y^^ , (43) 

the locus Z is empty and there is no emission. This critical speed depends on 
the surface tension of the fluid: for the case of a water/air interface it is equal to 
Vmin — 0.23 m/s. The absence of emission corresponds to a vanishing wave resis- 
tanee experienced by the slowly moving object which is able to slide with no fric- 
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tion on the surface of the fluid [10]. Still, the locaUzed deformation of the surface 
around the object is responsible for a renormalization of the mass of the object [9]. 

An efficient emission of surface waves with the associated wave resistance [10] is 
suddenly recovered as soon as v > Vmin. In this regime, the locus Z shown in Fig.7(e) 
has a kind of oval shape [9], with two intersections with the kx axis at respectively 



(1,2) 




(44) 



As expected, the two solutions merge to kx — ky for v > v^a, while at larger v » Vn 
they respectively tend to 



^1) 



2/„2 
min 



(45) 
(46) 



the former solution ki tends to zero in the large v limit and corresponds to almost 

(2) 9 

pure gravity waves, while the latter one kx quickly diverges as v and corresponds 
to almost pure capillary waves. 



6.2.1 Fast speed v » Vmin (deep fluid, negligible surface tension) 

Within the v ^ Vmin limit, we can start our discussion from the low wave vector 
region k <^ ky, where the waves have a mostly gravity nature, Q (k) ~ \/gk. Because 
of the fourth power of v/vmin that appears in (44), this limit is achieved already for 
moderate values of v/ Vmin of the order of a few unities. In this region, the locus E is 
approximately defined by the condition 



v^k^ 



K = ^[ki-°^)^ (47) 

whose shape is plotted in the panel (b) of Fig.7. The locus E extends in the \kx\ > 
k^ = gjv^ regions: for kx > ^, it has the form 



while for large kx ^ k^, one recovers an asymptotic behavior 



ky = ± — kl. (49) 
8 



A most remarkable feature are the inflection points at kf^ = ±i/3/2^/v^ where 
the slope dky/dkx is minimum. At these points the normal to the locus E makes 
the maximum angle to the kx axis, with a value (^max such that tan (^max = 1/V8~ 
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Fig. 7 Top row: dispersion of surface wave in the /i = oo deep water limit. Tlie dashed line indicates 
the i2 = k- V plane for generic particle speeds v = 2m/s (left (a) panel), v = 0.26m/s (right (d) 
panel). Middle row: coiTesponding shapes of the k-space locus E of resonantly excited modes 
(panels (b,e)); the green arrows indicate the normal to the locus Z, that is the direction of the relative 
group velocity v'^. Bottom row: real space patterns of the surface height modulation (panels (c,f)). 
These patterns are numerically obtained via a fast Fourier transfonn of the k-space perturbation 
(16) using the density and surface tension values of water The (c) panel corresponds to the Kelvin's 
ship-wave pattern behind a duck swimming on a deep lake, a picture of which is shown in the 
upper left panel of Fig.6. An original sketch by Lord Kelvin is shown in the lower panel of the 
same figure. 
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19°28'. This angle determines the aperture of the wake cone behind the moving 
object: remarkably, this value is universal and does not depend on the speed of the 
object. A picture of the Kelvin's ship-wave pattern behind a swimming duck on a 
quiet lake is shown in the upper left panel of Fig. 6; an original hand drawing by 
Lord Kelvin illustrating this physics is reproduced in the lower panel of the same 
figure. 

Another, related feature that is worth noticing is that for each angle |<^)| < <|>max, 
there exist two points on the locus E such that the normal to T, makes an angle 

with the negative kx axis: according to the theory discussed in Sec.2, these two 
solutions are responsible for the two interpenetrating fringe patterns: the so-called 
transverse waves with a small ky and the so-called diverging waves with large ky. 

The transverse waves are clearly visible as the long wavelength modulation right 
behind the object along the axis of motion: their wavevector fc;^ = fcj = ^/v^ is de- 
termined by the intersection of the locus T, with the kx axis. Remarkably, the faster 
the object is moving, the smaller is the wavevector kx- On Kelvin's hand drawing of 
Fig.6, the wavefront passing by point C belongs to the transverse wave pattern. 

The diverging waves are easily identified in the hand drawing as the wavefronts 
with opposite curvature connecting the source at A with the edge of the wake pattern 
where the two patterns collapse onto each other. The direction of the peculiar fringe 
modulation of the edge of the wake [indicated by the blue dashed lines on Fig. 7(c)] 
is determined by the wavevector k'"** of the inflection point of the k-space locus 
Z: the orientation of k'"" fixes the angle J3 to a value such that tan /3 = ky/kx\_^^f^ = 
I/V2, i.e. J3~35°. 

Of course, a complete treatment of the wake would require including the capil- 

(2) 

lary waves at very high wavevector k^v kx ^ ky, i.e. the part of the locus Z that 
closes the curve at large k's outside the field of view of Fig.7(b). However, the am- 
pUtude in these short-wavelength modes is quickly damped by viscous effects, so 
their contribution to the observable pattern turns out to be irrelevant in most practical 
cases. 



6.2.2 Moderate speed v > Vmin (deep fluid, significant surface tension) 

For moderate speeds v > Vmin, surface tension effects are no longer negligible and 
all points of the k-space locus Z contribute to the real-space pattern. In particular, 
the locus Z is a closed curve that does not encircle the origin, as shown in Fig. 7(e): 

(1 2) 

the two intersections with the kx axis at kx — kx , corresponding to gravity and 
capillary waves propagate with relative group velocities directed in opposite direc- 
tions from the fishing line of the celebrated experiment by Thomson. The pattern of 
long-wavelength gravity waves is located downstream of the fishing fine, while the 
short-wavelength capillary waves are located in the upstream region, see Fig.7(f) 
and the drawing by J. S. Russell reproduced in the middle panel of Fig.6. 

An approximate analytical understanding of this pattern can be obtained by ap- 
proximating the locus Z of Fig. 7(e) with a pair of circles analogously to the case of a 
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parabolic dispersion discussed in Sec.5 and shown in Fig.5(g-i): within this approx- 
imation, the shape of the wavefronts consists a system of hyperbolas, with a closer 
spacing ahead of the object. The qualitative agreement of the hyperbohc wavefronts 
of Fig.5(i) with the fuU calculations shown in Fig.7(f) is manifest. 

6.2.3 Effect of the source structure factor 

To complete the discussion, it is worth mentioning that the emission of waves can be 
hindered by the source structure factor 5(k) even at large v > Vmin- For example, the 
emission of surface waves will be strongly suppressed if the size £ of the source term 
(modelled as a Gaussian-shaped potential 5(r)) is large enough to have H » 1 for 
all points on L. For instance, for an object of typical size ^ = 30 cm, the argument 
such that ki^'^i < 1 imposes a lower critical speed v > vf'^ = 3m/s to the emission 
of gravity waves. A similar argument for capillary waves was mentioned to explain 
the characteristic swimming speed of some floating insects [49]. 



6.3 Shallow fluid 

When the wavelength of the perturbation is longer than the depth h of the fluid, the 
tanh(fc/z) term in the dispersion begins to be important and causes a radical change 
in the structure of the locus Z. The left and right columns of Fig. 8 illustrate the two 
regimes h > V^iy and h < \/3£y where the first correction to the sonic dispersion 
has either a sub-Unear or a super-linear nature. 

6.3.1 Small surface tension (sub-luminal dispersion) 

We start here from the case where the surface tension is small enough to have h > 
\/l)ly. Depending on the speed v of the object, several regimes can be identified. 
For very low speeds v < Vniin, the locus 1, is empty and there is no perturbation 
to the fluid. For intermediate speeds Vmin < v <C ^/^, the shape of the locus Z is 
determined by the high-k capillary region of the dispersion and is almost unaffected 
by the finite height h of the fluid: as in the deep fluid limit, 1, consists of closed, 
egg-shaped smooth curve and the real-space pattern again resembles a system of 
hyperbolas extending to the whole space, as illustrated in Fig.7(d-f). For increasing, 
yet sub-sonic speeds v < c., — -/gh, gravity recovers an important role, while the 
finite height h keeps providing only a small correction to the deep water behavior 
illustrated in Fig.7(a-c). 

The situation is completely different for supersonic speeds v > q [Fig.8(a-c)]: in 
this case, the locus Z starts from the origin, where it exhibits a conical singularity 
as a result of the sonic dispersion. At larger k, the locus T, recovers a shape similar 
to the infinitely deep fluid case: the sub-linear growth of the dispersion with k is 
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Fig. 8 Top row: dispersion of surface wave on shallow water of height h = 0.2m (left (a) panel) 
and h = 0.001 m (right (d) panel). The dashed line indicates the 12 = k- v plane for generic particle 
speeds v = 2m/s (left (a) panel), v = 0.14m/s (right (d) panel). Middle row: corresponding shapes 
of the k-space locus Z of resonantly excited modes (panels (b,e)). The dashed lines indicate the 
Cerenkov cone in the low wavevector region k^<^\. Bottom row: real space patterns of the surface 
height modulation (panels (c,f)). The black dashed lines indicate the Mach cone. These patterns are 
numerically obtained via a fast Fourier transform of the k-space perturbation (16) using the density 
and surface tension values of water The left panels correspond to a case where h > \/3(y and the 
lowest-order correction to the sonic dispersion (42) is sub-linear. The right panels correspond to 
a case where h < \/3ly and the lowest-order correction has the same super-linear behavior as the 
Bogoliubov dispersion (22) illustrated in Fig.4. 
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responsible for the fast increase of ky as a function of kx. Of course, the super- 
sonic dispersion of capillary waves at very large k [well outside the field of view 
of Fig.8(a,b)] makes the locus I, to close on itself. However, as already mentioned, 
these short-wavelength waves are quickly attenuated and hardly visible. 

The singularity of Z at the origin k = is responsible for the Mach cone and 
the disappearance of the transverse wave pattern, as shown in Fig. 8(c). As usual 
for sonic dispersions, the aperture of the Mach cone [indicated by the dashed line 
on Fig. 8(c)] depends on the source speed v as sin^ = Cj/v: on the T, locus shown 
in Fig. 8(b), this corresponds to the fact that the normal to L starts at a finite angle 
^ with the negative kx axis for k = 0. For growing k's, the angle monotonically 
decreases to meaning that the perturbation is restricted to the spatial region inside 
the Mach cone. The absence of the transverse wave pattern is clearly visible in 
Fig. 8(c) as the absence of modulation along the negative a: direction right behind the 
object. 

6.3.2 Shallow one-dimensional channel 

The restriction of this model to a one-dimensional geometry provides interesting 
insight on the physics of long wavelength surface waves in a spatially narrow chan- 
nel of width W. Spatial confinement along the orthogonal direction (say y) makes 
the corresponding wavevector to be quantized in discrete values determined by 
the boundary conditions at the edges of the channel, ky = np/W with the integer 
p = 0,1,2,. . .. For simpUcity, we assume that the transverse shape of the source 
(e.g. a boat saiUng along the channel) is broad enough to only excite the lowest 
mode at ky = 0, corresponding to a transversally homogeneous wave. 

Neglecting for simplicity also capillarity effects, a generahzed Landau criterion 
for one-dimensional gravity waves anticipates that a uniformly moving object can 
emit ky = gravity waves only if its velocity is slower than a maximum velocity 



This feature is easily understood looking at the ky = cut of the dispersion shown in 
Fig. 8(a): for sub-sonic speed v < c„ the ^2 = kxV straight Une corresponding to the 
Cerenkov condition (9) has a non-trivial intersection with the dispersion law, while 
the intersection reduces to the irrelevant kx = point for supersonic speeds v > Cj. 
In this case, no modes are any longer available for the emission, which reflects 
into a marked decrease of the wave drag friction experienced by the object. The 
observation of an effect of this kind when a ship travels at sufficiently fast speed 
long a channel was first reported by Scott Russel and often goes under the name of 
Houston paradox in the hydrodynamics and naval engineering literature [4]. 

As compared to the standard Landau criterion for superfluidity, it is interesting 
to note that the condition on the object speed to have a (quasi-)frictionless flow is 
here reversed: friction is large at slow speeds and suddenly drops for v > c^. This 



Vmax = max 




(50) 
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remarkable difference is due to the different sub- sonic rather than super- sonic shape 
of the gravity wave dispersion with respect to the Bogoliubov one. Of course, this 
suppression of friction is less dramatic when also higher p > transverse modes 
of the channel can be excited and a richer phenomenogy can be observed [50]. As 
we have previously discussed at length, in a transversally unlimited geometry the 
transition from sub-sonic to super-sonic speeds manifests itself as the disappearance 
of the transverse wave pattern from Kelvin's wake and a corresponding sudden but 
only partial decrease of the friction force. 

6.3.3 Large surface tension (super-luminal dispersion) 

In the opposite regime of large surface tension h < \fliiy, the super-linear form of 
the dispersion (42) makes the physics to closely resemble the behavior of impurities 
in a dilute superfluid discussed in Sec.4. For a slowly moving object at v < q, the 
locus 1, is empty and the perturbation of the surface remains localized in the vicinity 
of the impurity. For a fast moving object at v > Cj, the locus T, and the wake pattern 
closely resemble the corresponding ones for the case of a supersonically moving 
impurity in a superfluid shown in Fig.4: a Mach cone of aperture sin ^ = Cslv located 
behind the object [indicated by the dashed line on Fig. 8(f)] and a series of curved 
wavefronts ahead of the object. The most significant difference with the h > \/l>ly 
case of Sec. 6. 3.1 is the position of the modulation with respect to the Mach cone: in 
the sub-linear case of panel (c), it lies within (i.e. behind) the Mach cone, while in 
the super-linear case of panel (f), it stays outside (i.e. in front of) the Mach cone. 



7 Cerenkov processes and the stability of analog black/white 
holes 

The systems that were considered in the previous sections are presently among the 
most promising candidates for the realization of condensed matter analogs of grav- 
itational black (or white) holes: the key idea of analog models is to tailor the spatial 
structure of the flow in a way to obtain a horizon surface that waves can cross only in 
one direction. Upon quantization, a number of theoretical works have predicted that 
a condensed matter analog of Hawking radiation should be emitted by the horizon. 
A complete review of this fascinating physics can be found in the other chapters 
of the book. In this last section, we shall review some consequences of Cerenkov 
processes that are most significant for the stabiUty of analog black and white holes 
based on either flowing superfluids or surface waves on flowing water. 

The role of Cerenkov-like emission processes in the dynamics of the strong op- 
tical pulses that are used in optical analog models based on nonUnear optics [5 1] 
is still in the course of being elucidated and interesting experimental observations 
in this direction have recently appeared [52]. Here, it is important to remind that, 
differently from the all-optical Cerenkov radiation experiments of [17, 31] where 
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an effective moving dipole was generated by nonlinearity, the analog models 
of [51] are based on the time- and space-dependent effective refractive index profile 
due to a x^^^ optical nonlinearity: given the centro-symmetric nature of the medium 
under examination, no effective moving dipole can in fact appear unless the medium 
shows some material imperfection. 



7.1 Superfluid-based analog models 

Let us start from the simplest case of analog blaclc/white holes based on flowing 
superfluids for which a complete theoretical understanding is available [25, 53, 54]. 
In a one dimensional geometry, the horizon consists of a point separating a region 
of sub-sonic flow from a region of super-sonic flow. In a black hole the sub-sonic 
region hes upstream of the horizon, while in a white hole the sub-sonic region hes 
downstream of the horizon. A sketch of both configurations is reproduced in Fig.9 
together with the Bogoliubov dispersion of excitations as observed from the labora- 
tory frame: in the most common configurations, the flow has a non-trivial structure 
only in a small region around the horizon and recovers a homogeneous shape with 
space-independent density and speed farther away from the horizon. In the labora- 
tory reference frame (corresponding to the rest frame of the impurity), the Cerenkov 
emission occurs in the zero-frequency BogoUubov modes, the so-called zero modes. 

As we have discussed in detail in Sec. 4, the super-linear nature of Bogoliubov 
dispersion restricts Cerenkov emission processes to super-sonic flows, where they 
generate waves that propagate in the upstream direction. In the geometry under 
consideration here, the flow is everywhere smooth exception made for the hori- 
zon region. Combining these requirements immediately rules out the possibility of 
Cerenkov emission in black hole configurations: the group velocity of the zero mode 
waves emitted at the horizon points in the direction of the sub-sonic region, where 
it can no longer be supported. This simple kinematic argument contributes to ex- 
plaining the remarkable dynamical stabiUty of acoustic black hole configurations, 
as observed in numerical simulations of their formation starting from a uniformly 
moving fluid hitting a localized potential barrier [55]. 

In contrast, white hole configurations are much more sensitive to the dissipation 
of energy via Cerenkov processes: Bogoliubov excitations can appear in the super- 
sonic region upstream of the horizon and give rise to significant modulations of 
the density and flow speed, the so-called undulation patterns. Several reasons make 
such processes to be potentially harmful to the study of quantum features of the 
white hole radiation. To the best of our knowledge, the only known realistic scheme 
to generate a white hole configuration in a flowing superfluid is the one of [25] using 
a simultaneous spatial and temporal modulation of both the atom-atom interaction 
strength and the external confining potential. The main difficulty of this configura- 
tion is that it requires a very precise tuning of the system parameters to eUminate 
unwanted Cerenkov emission processes that may mask the quantum vacuum radia- 
tion. 
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Fig. 9 Main panels: sketch of the flow geometry for white (a) and black (b) hole contigurations 
based on superfluids. Smaller panels (al,a2,bl,b2): dispersion of Bogoliubov excitations in the 
asymptotic regions far from the horizon. Cerenkov emission is only possible for white hole con- 
figurations: the con'esponding zero mode is indicated in blue in (al). Figure from [25]. 
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Fig. 10 Main panels: sketch of the flow geometry for white (a) and black (b) hole configurations 
based on surface waves on a fluid. Smaller panels (al,a2,bl,b2): surface wave dispersion in the 
asymptotic regions far from the horizon. Cerenkov emission is only possible for white hole config- 
urations: the corresponding zero mode is indicated in blue in (a2). Differently from the superfluid 
case of Fig. 9, the zero mode now has a group velocity in the downstream direction into the sub- 
sonic region. Parameters: flow speed v = 2 m/s, fluid depth /i = 0.1 m [white hole, upstream inner 
region, panel (al)], v = 0.666 m/s, fluid depth h = 0.3 m [white hole, downstream outer region, 
panel (a2)]; flow speed v = —2 m/s, fluid depth /i = 0.1 m [black hole, downstream inner region, 
panel (bl)], v = —0.666 m/s, fluid depth h = 0.3 m [black hole, upstream outer region, panel (b2)]. 
For the chosen parameters, the effect of surface tension at the water/air interface is negligible. 
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Even if a perfect preparation of the white hole is assumed, Cerenkov emission 
processes may still be triggered by nonlinear effects in the horizon region. As it 
was shown in [25], an incident classical Bogohubov wavepacket is able to induce a 
distortion of the horizon proportional to the square of its amplitude, which then re- 
sults in the onset of a continuous wave Cerenkov emission from the horizon and the 
appearance of a spatially oscillating modulation in the density profile upstream of 
the horizon. Of course, a similar mechanism is expected to be initiated by quantum 
fluctuations when back-reaction effects are included in the model, i.e. the non-linear 
interaction of quantum fluctuations with the underlying flow. A third, more subtle 
mechanism of instability of a white hole configuration was unveiled in [25]: the 
1 divergence of the matrix elements of the S-matrix for low-frequency outgo- 
ing modes in the neighborhood of the finite wavevector zero mode is responsible for 
a steady growth of the density fluctuation ampUtude in time since the formation of 
the white hole. Even if the temporal growth of fluctuations follows a slow logarith- 
mic (hnear at a finite initial temperature J > 0) law, still it is expected to strongly 
affect the properties of the horizon at long times. 

The situation is expected to be different if fully three dimensional systems with- 
out transverse confinement are considered. In this case, Cerenkov emission can take 
place also in a black hole configuration: a distortion of the horizon by classical or 
quantum fluctuations with a non-trivial transverse structure is in fact able to excite 
Bogoliubov modes with a finite transverse component ky ^ 0. As we have seen in 
Fig.4, there exist such modes that can propagate in the downstream direction into 
the supersonic region inside the black hole. The effect of these zero modes on the 
finite-^y S-matrix for the black hole is presently under investigation [56]. 



7.2 Analog models based on surface waves 

The different dispersion of surface waves is responsible for dramatic qualitative 
differences in the wave propagation from the horizon of analog black and white 
holes configurations. As it is sketched in Fig.lO(a,b), the trans-sonic interface is 
generally created in these systems by means of a spatial variation of the fluid depth 
h [57, 58]. Depending on the detailed shape of the transition region and/or on the 
presence of fluctuations, Cerenkov emission by the horizon can occur into the zero 
modes of zero energy in the laboratory frame, as observed e.g. in [59]: a quantitative 
estimate of the amplitude of the resulting stationary undulation pattern for specific 
configurations of actual experimental interest requires however a complete solution 
of the hydrodynamic equations, which goes beyond the scope of the present work. 

With an eye to the experiments of [57, 58], we can restrict our attention to 
low-wavevector gravity waves and neglect surface tension effects in a simplest 
one-dimensional geometry. In this case, Cerenkov emission processes only oc- 
cur for flow speeds lower than the sonic speed Cs = \/gh and result in an emis- 
sion in the downstream direction. From the surface wave dispersions shown in 
Fig.l0(al,a2,bl,b2), one can easily see that Cerenkov emission can again only oc- 
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cur in white hole configurations, which are then expected to be again less stable 
than black hole ones. As in the superfluid case, Cerenkov emission of surface waves 
with a finite ky^Q become possible also in the black hole case as soon as a two 
dimensional geometry is considered and Kelvin's diverging waves are allowed. In 
spite of its importance in view of on-going experiments, we are not aware of any 
comprehensive work having studied in full detail the dynamical stabihty of surface 
wave-based analog white/black holes as in the case of superfluid-based ones [25]. 



8 Conclusions 

In this chapter we have presented a review on some most significant aspects of the 
Cerenkov effect from a modern and interdisciplinary point of view. In our perspec- 
tive the same basic process of generalized Cerenkov emission encompasses all those 
emission processes that take place when a uniformly moving source is coupled to 
some excitation field: as soon as the source velocity exceeds the phase velocity of 
some mode of the field, this gets continuously excited. Simple geometrical argu- 
ments are presented that allow to extract the shape of the emission pattern in real 
and k space from the dispersion law ^2 (k) of the field. Application of the general 
concepts to some most illustrative cases is discussed, from the standard Cerenkov 
emission of relativistically moving charged particles in non-dispersive media, to the 
Mach cone behind a supersonically moving impurity in a superfluid, to the wake 
of gravity and capillary waves behind a duck swimming on the surface of a quiet 
lake. The impact of Cerenkov emission processes on condensed-matter analog mod- 
els of gravitational physics is finally discussed. Open questions in this direction are 
reviewed. 
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